Tuning the conductance of Dirac fermions on the surface of a topological insulator 
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We study the transport properties of the Dirac fermions with Fermi velocity vf on the surface of a 
topological insulator across a ferromagnetic strip providing an exchange field J over a region of width 
d. We show that the conductance of such a junction changes from oscillatory to a monotonically 
decreasing function of d beyond a critical J . This leads to the possible realization of a magnetic 
switch using these junctions. We also study the conductance of these Dirac fermions across a 
potential barrier of width d and potential Vo in the presence of such a ferromagnetic strip and show 
that beyond a critical J , the criteria of conductance maxima changes from x = eVod/hvp = rm to 
X = {n-\-l/2)n for integer n. We point out that these novel phenomena have no analogs in graphene 
and suggest experiments which can probe them. 

PACS numbers: 71.10.Pm, 73.20.-r 
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Topological insulators in both two- and three- dimen- 
sions (2D and 3D) have attracted a lot of theoretical and 
experimental attention in recent years P, |3, S iJ| • It has 
been shown in Rcf. [4] that such 3D insulators can be 
completely characterized by four integers z^o and vi,2,'i- 
The former specifies the class of topological insulators to 
be strong (i/q = 1) or weak [uq — 0), while the latter in- 
tegers characterize the time-reversal invariant momenta 
of the system given by Mq = {vibi, t^2^2, t^3^3)/2, where 
^1,2,3 are the reciprocal lattice vectors. The topological 
features of strong topological insulators (STI) are robust 
against the presence of time-reversal invariant perturba- 
tions such as disorder or lattice imperfections. It has been 
theoretically predicted and experimentally verified 

that the surface of a STI has an odd number of Dirac 
cones whose positions are determined by the projection 
of Afo on to the surface Brillouin zone. The position 
and number of these cones depend on both the nature of 
the surface concerned and the integers J^i,2,3- For several 
compounds such as HgTe and Bi2Se3, specific surfaces 
with a single Dirac cone near the F point of the 2D Bril- 
louin zone have been found [2, Q . Such a Dirac cone is 
described by the Hamiltonian 
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where f?(/) denotes the Pauli (identity) matrices in spin 
space, ij} = (V't, ^x)"^ is the annihilation operator for the 
Dirac spinor, vp is the Fermi velocity, and /i is the chemi- 
cal potential [6] . Recently, several novel features of these 
surface Dirac electrons such as the existence of Majorana 
fermions in the presence of a magnet-superconductor in- 
terface on the surface d, 0, generation of a time- 
reversal symmetric Px + ipy superconducting state via 
proximity to a s-wave superconductor j^], anomalous 
magnetoresistance of ferromagnet-ferromagnet junctions 
[i| and novel spin textures with chiral properties [13] have 
been studied in detail. 




FIG. 1: Proposed experimental setups: a) Left panel: The 
ferromagnetic film extends over region II of width d provid- 
ing an exchange field in this region, b) Right panel: The 
film extends over region III while the region II has a barrier 
characterized by a voltage Vb- V and J denote the bias volt- 
age and current across the junction respectively. See text for 
details. 



In this letter, we study the transport properties of 
these surface Dirac fermions in two experimentally re- 
alizable situations shown in Fig. [TJ The first study con- 
cerns their transport across a region with a width d with 
a proximity-induced exchange field arising from the 
magnetization rfi — moy of a proximate ferromagnetic 
film as shown in the left panel of Fig. [T] We demon- 
strate that the tunneling conductance G of these Dirac 
fermions through such a junction can either be an oscilla- 
tory or a monotonically decaying function of the junction 
width d. One can interpolate between these two quali- 
tatively different behaviors of G by changing toq (and 
thus J) by an applied in-plane magnetic field leading to 
the possible use of this junction as a magnetic switch. 
The second study concerns the transport properties of 
Dirac fermions across a barrier characterized by a width 
d and a potential Vq in region II with a magnetic film 
proximate to region III as shown in the right panel of 
Fig. [H We note that it is well known from the context 
of Dirac fermions in graphene Uj that such a junction. 
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in the absence of the induced magnetization, exhibits 
transmission resonances with maxima of transmission at 
X — eVod/Hvp — TiTT, where n is an integer. Here we 
show that beyond a critical strength of mo, the maxima 
of the transmission shifts to x = (rt + l/2)7r. Upon fur- 
ther increasing toqi one can reach a regime where the 
conductance across the junctions vanishes. We stress 
that the properties of Dirac fermions elucidated in both 
these studies are a consequence of their spinor structure 
in physical spin space, and thus have no analogs for ei- 
ther conventional Schrodinger electrons in 2D or Dirac 
electrons in graphene [l^ . 

We begin with an analysis of the junction shown in the 
left panel of Fig. [1] The Dirac fermions in region I and 
III are described by the Hamiltonian in Eq. ([T]). Con- 
sequently, the wave functions of these fermions moving 
along ±x in these regions for a fixed transverse momen- 
tum ky and energy e can be written as 
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where i takes values I and III, a = a,vcsh\{hv pky / \ e -f 
and kx{e) = y^[(e -I- ^)/hvFY — ky. In region II, the 
presence of the ferromagnetic strip with a magnetiza- 
tion mo = ruQy leads to the additional term iJinduccd = 
/ dxdyje{x)9{d - x)lP^x)ay ip{x), where mo is 

the exchange field due to the presence of the strip 
and 9{x) denotes the Heaviside step function. Note 
that -ffinducod may be thought as a vector potential 
term arising due to a fictitious magnetic field Bf = 
{J / evF)[5{x) — S{d — x)]z. This analogy shows that our 
choice of the in-pane magnetization along y is completely 
general; all gauge invariant quantities such as transmis- 
sion are independent of x-component of mo in the present 
geometry. For a given mo, the precise magnitude of J 
depends on several factors such as the exchange coupling 
of the film and can be tuned, for soft ferromagnetic films, 
by an externally applied field The wave function for 
the Dirac fermions in region II moving along ±x in the 
presence of such an exchange field is given by 



(3) 



where (3 — a.rcsm{hvF{ky + M)/\e + ^\), M — J j {hvp), 
and fc^(e) = ^/[(e + ^j)/hvFf - {ky + M)^. Note that 
beyond a critical Mc = ±2|e -I- ^\/{hvF) (and hence a 
critical Jc = ±2|e -I- k^ becomes imaginary for all ky 
leading to spatially decaying modes in region II. 

Let us now consider an electron incident on region II 
from the left with a transverse momentum ky and energy 
e. Taking into account reflection and transmission pro- 
cesses at X = and x — d, the wave function of the elec- 
tron can be written as ipi — tpf+rtpj , ip^^ = pip^j+qtpjj, 
and ipiii — tip'j^j]-. Here r and t are the reflection and 
transmission amplitudes and p (q) denotes the amplitude 
of right (left) moving electrons in region II. Matching 
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FIG. 2: Plot of tunneling conductance G/Go for a fixed V 
and ^ as a function of the effective width z = d\eV + fi\/hvF 
for hvFM/\eV + mI = 0.3 (green dotted line), 0.7 (blue solid 
line), 1.3 (black dash-dotted line) and 2.1 (red dashed line). 
The value of the critical AI is given by hvFM/\eV + nl = 2. 
See text for details. 



boundary conditions on ijjj and ipji at x = and ipn 
and tpiii at X = d leads to 



1 + r — p 



te 



^^i{k^d+a) ^ pgi{k^d+l3) _ ^^-i{k^d+l3) ^ 



(4) 



Solving for t from Eq. ^ , one finally obtains the con- 
ductance G = dl/dV = (Go/2) /^//2 ^ cos(a)da. Here 

Go = p{eV)we'/{iTh^VF), p{eV) ^]{tJi + eV)\/[2n{hvFf] 
is the density of states (DOS) of the Dirac fermions and 
is a constant for p, 3> eV , w is the sample width, and the 
transmission T — |ip is given by 

T = cos^(q;) cos^(/3)/[cos^(/c^(i) cos^(a) cos^(/3) 

+ sin2(/c^d)(l -sin(a)sin(/3))2]. (5) 

Eq. ([5]) and the expression for G represent one of the 
main results of this work. We note that for a given a, T 
has an oscillatory (monotonically decaying) dependence 
on d provided k^ is real (imaginary). Since k^ depends, 
for a given a, on M , we find that one can switch from an 
oscillatory to a monotonically decaying d dependence of 
transmission in a given channel (labeled by ky or equiva- 
lently a) by turning on a magnetic field which controls mo 
and hence M. Also since — 1 < sin(Q;) < 1, we find that 
beyond a critical M — Mc, the transmission in all chan- 
nels exhibits a monotonically decaying dependence on d. 
Consequently, for a thick enough junction one can tune G 
at fixed V and from a finite value to nearly zero by tun- 
ing M {i.e., mo) through Mc. Thus such a junction may 
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FIG. 3: Plot G/Go versus eV/fi for several representatives 
values hvpM/fi ranging from 3 (left-most black solid curve) 
to 5 (right-most magenta dash-double dotted line) in steps of 
0.5 The effective junction width zq = 5 for all plots. The inset 
shows a plot of eVc/fJ^ versus hvpM/fi. See text for details. 



be used as a magnetic switch. These qualitatively differ- 
ent behaviors of the junction conductance G for M below 
and above Mc is demonstrated in Fig. [2] by plotting G as 
a function of effective barrier width z — d\eV + fj,\/hvF 
for several representative values of hvpM /\eV + fi\. Since 
T and hence G depends on M through the dimcnsionless 
parameter hvpM/\eV + this effect can also be ob- 
served by varying the applied voltage V for a fixed /i, d, 
and M. In that case, for a reasonably large dimcnsionless 
barrier thickness zq = dfj,/hvF, G/Gq becomes finite only 
beyond a critical voltage \eVc + fJ.\ — hvFM/2 as shown 
in Fig. [3] for several representative values of zq. This crit- 
ical voltage Vc can be determined numerically by finding 
the lowest voltage for which G/Gq exhibits a monotonic 
decay as a function of zq. The plot of eVc/ fJ- as a function 
of hvpM/fi, shown in inset of Fig. [31 demonstrates the 
expected linear relationship between Vc and M. We note 
such a magnetic field or applied bias voltage dependence 
of the junction conductance necessitates that the Dirac 
electrons represents spinors in physical spin space and is 
therefore impossible to achieve in graphene . 

Next, we analyze the junction shown in the right panel 
of Fig. [T] where the region III below a ferromagnetic film 
is separated from region I by a potential barrier in re- 
gion II. Such a barrier can be applied by changing the 
chemical region in region II either by a gate voltage Vq 
or via doping fH*]. In the rest of this work, we will analyze 
the problem in the thin barrier limit for which Vq —^ oo 
and d ^ 0, keeping the dimensionless barrier strength 
X = eVod/{hvp) finite. The wave function of the Dirac 
fermions moving along ±x with a fixed momentum ky 




FIG. 4: Plot of tunneling conductance Gi/Go versus the ef- 
fective barrier strength x ^-nd hvpM for fixed applied voltage 
V and chemical potential ^j,. Gi vanishes for \M\ > Mc = 
2\eV + n\/hvF. 

and energy e in this region is given by 

ij'jf = (l,±e=^*^) e'(±'=^'^+''«^V^/2, (6) 

where 7 = arcsin^hv pky/\e + eVo + and k^{e) = 
y^[(e + cVq + /i) /hvp]'^ — fc^. The wave functions in re- 
gion I and III are given by Eqs. ^ and ^ respectively: 
"0/ = V'/ snd V'/// — i^ii- Note that one can have a 
propagating solution in region III only if \M\ < \Mc\- 

The transmission problem for such a junction can be 
solved by an procedure similar to the one outlined above 
for the magnetic strip problem. For an electron ap- 
proaching the barrier region from the left, we write down 
forms of the wave function in the three regions I, II 
and III: ip\ ^ + n^Y , tp'jj = Pii^jf + gW^/f, and 
= ti'ipu- outlined earlier, one can then match 
boundary conditions at x — and x = d, and obtain the 
transmission coefficient Ti — \ti\'^k^/kx as 

Ti = 2cos(/3)cos(a)/[l +cos(/3- a) 

- cos^(x){cos(/3 - a) - cos(/3 + a)}]. (7) 

Note that in the absence of the ferromagnetic film over 
region III, /3 = a, and Ti ^ cos2(a)/[l - 

cos (x)sin2(a)]. The expression for T^, reproduced 
here for the special case of A/ = 0, is well known 
from analogous studies in the context of graphene, and 
it exhibits both Klein paradox (T" — 1 for a — Q) 
and transmission resonance {T^ = 1 for x = 
[13]. When M 7^ 0, we find that the transmission 
for normal incidence {ky = 0) does become indepen- 
dent of the barrier strength, but its magnitude deviates 
from unity: Tf°™'^i = 2^1 - {hv p M / \eV + ^l\Y / {I + 
y/1 ~ {hvpM/\eV + The value of rn°"°ai de- 

creases monotonically from 1 for M = to for | = 
\eV + fi\/{hvp) and can thus be tuned by changing M 
(or V) for a fixed V (or M) and /i. 
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FIG. 5: Top left panel: Plot of Gi/Go versus x for 
hvFM/\eV + nl = 0.1 (black solid line), 0.7075 (red dashed 
line), and 1.4 (blue dash-dotted line) for fixed V and /i. Top 
right panel: Plot of Gi/Go versus x a,t M = M* showing 
the period halving. Bottom left panel: Plot of Xmax versus 
hvpM /\eV + /^l showing the shift of conductance maxima. 
Bottom right panel: Plot of AGi/Go versus hvpM /\eV + /i| 
which crosses for M < Mc at M = M*. The dotted line is 
a guide to the eye. 



The conductance of such a junction is given by Gi = 
(Go/2) /"^ Ti cos(a)da, where ai 2 are determined from 
the solution of cos(/3) = for a given M (Eq. ([3])). A plot 
of Gi as a function of hvpM /\eV + /i| and x (for a fixed 
eV and /i) is shown in Fig. 21 We find that the amplitude 
of Gi decreases monotonically as a function of \M\ reach- 
ing at M = Mc beyond which there are no propagating 
modes in region III. Also, as we increase M, the conduc- 
tance maxima shifts from x — ^''^ to x = (f^ + l/2)7r 
beyond a fixed value of M*{V) ~ zbcoleV + iM\/{hvF) as 
shown in top left panel of Fig. [S] Numerically, we find 
Co = 0.7075. At M = M*, Gi(x = mr) = Gi(x = 
[n -\- l/2)7r), leading to a period halving of Gi(x) from tt 
to 7r/2 . This is shown in top right panel of Fig. [5] where 
Gi(M — M*) is plotted as a function of x- We note that 
near M* , the amplitude of oscillation of Gi as a function 
of X becomes very small so that Gi is almost indepen- 
dent of X- In the bottom left panel of Fig. [5l we plot 
X = Xmax (the value of x at which the first conductance 
maxima occurs) as a function of hvpM /\eV + /i| which 
clearly demonstrates the shift. This is further highlighted 
by plotting AGi — Gi(x = 0) — Gi(x = 7r/2) as a func- 
tion oi hvpM l\eV + ^\ in the bottom right panel of Fig. 
[51 For M < Mc, AGi crosses zero at M = M* indicat- 
ing the position of the above-mentioned period halving. 
Thus we conclude that the position of the conductance 
maxima depends crucially on hvpM/\eV + /i| and can be 
tuned by changing either AI or V. 

The experimental verification of our results would in- 
volve preparation of junctions by depositing ferromag- 
netic films on the surface of a topological insulator. For 



the geometry shown in the left panel of Fig.[l] we propose 
measurement of G as a function of mo whose magnitude 
and direction can be tuned by an externally applied in- 
plane magnetic field for soft ferromagnetic films We 
predict that depending on mo, G should demonstrate ei- 
ther a monotonically decreasing or an oscillatory behav- 
ior as a function of d. Another, probably more experi- 
mentally convenient, way to realize this effect would be 
to measure Vc of a junction of width d for several M and 
confirm that Vc varies linearly with M with a slope of 
hvp/{2e), provided /x and d remain fixed. For the geom- 
etry depicted in the right panel of Fig. [U one would, in 
addition, need to create a barrier by tuning the chemical 
potential of an intermediate thin region of the sample as 
done earlier for graphene 1J| . Here we propose measure- 
ment of Gi as a function of Vq (or equivalently x) for 
several representative values of mo and a fixed V. We 
predict that the maxima of the tunneling conductance 
would shift from x — ni: to x = {n + l/2)7r beyond a 
critical mo for a fixed V , or equivalently, below a critical 
V ^ for a fixed mo. 

In conclusion, we have studied the transport of Dirac 
fermions on the surface of a topological insulator in the 
presence of proximate ferromagnetic films in two experi- 
mentally realizable geometries. Our study unravels novel 
features of the junction conductances which have no ana- 
log in either graphene or 2D Schrodinger electrons and 
can be verified in realistic experimental setups. 
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